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We have searched for proton decays via p — > e + 7r° and p — > fi + ir° using data from a 
91.7 kiloton-year exposure of Super-Kamiokande-I and a 49.2 kiloton-year exposure of Super- 
Kamiokande-II. No candidate events were observed with expected backgrounds induced by atmo- 
spheric neutrinos of 0.3 events for each decay mode. From these results, we set lower limits on the 
partial lifetime of 8.2 x 10 33 and 6.6 x 10 33 years at 90% confidence level for p — » e + 7r° and p — > /i + 7r° 
modes, respectively. 
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One of the unique predictions of Grand Unified The- 
ories (GUTs) is baryon number violation. Via the ex- 
change of a very heavy gauge boson, two quarks in a pro- 
ton can transform into a lepton and an anti-quark result- 
ing in a lepton plus meson final state. Because the mass 
of the gauge boson is predicted to be very heavy such 
processes are expected to be very rare. The prototypical 
GUT, minimal SU(5)0,HI|], predicts the lifetime of the 
proton when it decays to e + 7r° to be less than 10 32 years. 
A number of other GUT models favor proton decay into 
this final state 0, [H, or with comparable branching into 
the similar final state /^ + 7r o 0, Contemporary theo- 
ries must take steps to evade the stringent limits set by 
Super-Kamiokande [H and prior experiments^, |T(| that 
decisively ruled out minimal SU(5). One means is by in- 
troducing supersymmetry, which raises the mass of the 
heavy gauge boson, but which introduces other decay 
channels which are also constrained by experiment & 

In this paper, we describe our search for proton de- 
cay by the reactions p — ► e + 7r° and p — > p, + Tr°. The 
result of the search was negative, no proton decay was 
found. We update our lifetime limit on the e + 7r° channel 
based on six times greater exposure than our previous 
publication Q; we report for the first time our limit on 
the p, + ir° channel. Approximately 1/3 of the exposure 
was with photocollection coverage reduced by a factor of 
two, so we also describe the effect of reduced photocov- 
erage on these studies. 

Super-Kamiokande is a 50-kiloton water Cherenkov de- 
tector located in the Kamioka Observatory in Japan [l^]. 
The 22.5-kiloton fiducial volume of the detector was 
viewed by 11146 20- inch diameter photomultiplier tubes 
(PMTs) . Super-Kamiokande started observation in April 
1996 and stopped in September 2001 (SK-I) for a detec- 
tor upgrade. On November 12, 2001, during water filling 
after the upgrade, a shock wave was initiated by an im- 
ploding bottom PMT. About half of the PMTs were lost 
in this accident. In the first recovery stage (SK-II), 5182 
PMTs were enclosed in a fiber reinforced plastic case with 
an acrylic cover for protection against another chain re- 
action. SK-II re-started observation in October 2002 and 
stopped in October 2005 for full detector reconstruction. 
The transparency and the reflection of the acrylic cover 
were 97% and 1%, respectively. The photo-sensor cover- 
age was 40% and 19% in SK-I and SK-II, respectively. 

A Monte Carlo simulation was used to estimate the 
efficiency of detecting proton decay occurring in water 
(H2O). The two free protons and eight bound protons in 
an H 2 molecule are assumed to decay with equal prob- 
ability. For the case of a free proton in hydrogen, the 
momenta of the decay particles are uniquely determined 
by two-body kinematics. For the case of a bound proton 
in oxygen, the decay particle momenta are no longer de- 
termined by simple two-body decay; this is due to Fermi 
motion of the protons and the nuclear binding energy, as 
well as meson-nuclear interactions in oxygen. 



The Fermi momentum distribution of the nucleons was 
taken to be the same as that measured by electron scat- 
tering on 12 C for S-states and P-states[13]. Nuclear bind- 
ing energy was taken into account by using a modified 
proton mass given by M' p = M p — Eb, where M p is 
the proton rest mass, and Eb is the nuclear binding en- 
ergy. The value of Eb for each simulated event was ran- 
domly selected from a Gaussian probability density func- 
tion with O, cr) = (39.0, 10.2) MeV for the S-state and 
(fi, a) = (15.5, 3.82) MeV for the P-state. Pions inter- 
act strongly in the nucleus and may undergo scattering, 
charge exchange, or absorption as they travel through the 
nucleus. These interactions affect our ability to recon- 
struct the p — > e + 7r° event, and were carefully simulated 
to give our best estimate of the detection efficiency of our 
expected signal. The cross sections used in the simulation 
of each type of pion-nucleon interaction were calculated 
by the model of Oset et al. [l4[ • Another effect in the nu- 
clear medium which we take to occur for approximately 
10% of proton decay events is that of correlated decay: 
pN — > e + 7r°A [15]. The value of the invariant mass cal- 
culated in the case of correlated decays is lower than that 
of standard (uncorrelated) proton decays, leading to a re- 
duction in the proton decay detection efficiency when an 
invariant mass selection criterion is applied. 

Backgrounds to the proton decay search arise from at- 
mospheric neutrino interactions which may mimic the 
p — > e + 7r° and p — > /i + 7r° signals by charged current 
interactions such as vN — > £N'ir° and neutral current 
interactions such as vN vN'tt(tt's). Backgrounds 
were estimated using the neutfljj neutrino interaction 
Monte Carlo simulation with an input atmospheric neu- 
trino flux 17 1 . 

Particles produced in the simulation of proton decay 
events and atmospheric neutrino interactions were passed 
through a GEANT-3 [l8| based custom detector simu- 
lation to model Cherenkov light emission from charged 
particles, particle and light propagation through matter, 
detector geometry, and the response of the PMTs and 
data acquisition electronics. Hadronic interactions were 
treated by CALOR [l9| for nucleons and charged pions of 
Pjt> 500 MeV/c, and by a custom simulation program 
[201 ] for charged pions of p^ < 500 MeV/c. 

We used data from a 91.7 kiloton-year exposure of 1489 
live days during SK-I and a 49.2 kiloton-year exposure of 
798 live days during SK-II. There were about 10 6 event 
triggers/day above a threshold of a few MeV. We applied 
several stages of data reduction in order to remove back- 
ground from cosmic rays, flashing PMTs and radioactiv- 
ities. In SK-I and SK-II, 12232 and 6584 events were 
obtained, respectively. The background contamination 
of non-neutrino interactions in the final fiducial volume 
sample was negligible, less than 1%. The efficiency for 
proton decay events to appear in this sample was esti- 
mated to be greater than 99%. 

For events after the chain of reduction processes, physi- 
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cal quantities were reconstructed using timing and charge 
information from each PMT. Software reconstruction al- 
gorithms are shared between SK-I and SK-II analyses, 
taking into account the different PMT densities. The di- 
rections and opening angles of visible Cherenkov rings 
were reconstructed with a likelihood method, using the 
fitted vertex position where the time-of-flight subtracted 
timing distribution was most sharply peaked. Each iden- 
tified Cherenkov ring was categorized as a showering 
electron-like (e-like) or a non-showering muon-like (fi- 
like) ring. Momentum was determined by the sum of 
photoelectrons corrected for light attenuation in water, 
PMT angular acceptance, and PMT coverage. Finally, 
the number of electrons from muon decay was determined 
by searching for delayed electron signals. The delayed 
signals were found as delayed timing peaks within the 
primary event (up to 800 ns from the trigger time) , or as 
subsequent event triggers within 20 /js after the primary 
event. In multi-ring events such as proton decay signals, 
each PMT observes an integrated sum of Cherenkov light 
coming from all rings. In order to determine the momen- 
tum for each ring, the fraction of the photoelectrons from 
each ring was estimated for each PMT using an expected 
Cherenkov light distribution for each ring. Those pho- 
toelectrons were corrected by the amount of scattered 
Cherenkov light in water, reflected light on PMT sur- 
faces, and the vertex position shift due to the 7's conver- 
sion length. 

For SK-I (SK-II) free proton decays of p — > e + 7r°, the 
estimated vertex resolution was 18.1 (20.1) cm, the frac- 
tion of 2-ring and 3-ring events was 39±2% and 60±2% 
(38±2% and 60±2%), and the [i/e particle misidenti- 
fication probability was estimated to be 3.3% (3.4%). 
We define the total momentum in an event as Ptot — 

Eall rinqs -> i i -> . , , n , 

, Pi\, where pi is reconstructed momentum vec- 

tor of the i-th ring. The total invariant mass is also 
defined as M to t — \J Ef ot — P 2 ot , where total energy 

Etot = J2i nn9S \Z\Pi\ 2 + m i > an d m i i s assumed to be 
the electron (muon) rest mass for e-like (/x-like) ring. The 
total invariant mass for simulated free proton decays was 
well reconstructed both in SK-I and SK-II. The mean 
of reconstructed masses agreed within 1% between SK-I 
and SK-II. The resolution of the total momentum and 
total invariant mass was 30.5 MeV/c and 28.7 MeV/c 2 
for SK-I, and 36.6 MeV/c and 38.4 MeV/c 2 for SK-II, 
respectively. 

The energy scale stability of the detector was tested 
by the mean reconstructed energy of stopping cosmic 
ray muons and their decay electrons. It varied within 
±0.88% (±0.55%) over the SK-I (SK-II) period. The ab- 
solute energy scales for SK-I and SK-II were separately 
adjusted by using observed photo electrons of penetrating 
cosmic ray muons in each data taking period. The scale 
was checked by many calibration data such as stopping 
cosmic ray muons, electrons from their decays, and the 
invariant mass of 7r°s produced in atmospheric neutrino 



TABLE I: Number of surviving events in the data and the at- 
mospheric neutrino MC at each selection criterion in 91.7(SK- 
I)+49.2(SK-II) kiloton-year exposure. Detection efficiencies 
using the proton decay MCs are also listed. 



Criteria 


p -> e+n° 


p -> h + ty° 


Efficiency(%) 




data (atm.MC) 


data (atm.MC) 


(e+Tr ) 


(» + *°) 


in fiducial 


18816 (19269) 


18816 (19269) 


98.6 


98.8 


(A) 


4889 (5124) 


4889 (5124) 


73.1 


73.6 


(B) 


3036 (3141) 


1536 (1604) 


64.7 


61.4 


(C) 


2541 (2613) 


1281 (1284) 


62.6 


59.7 


(D) 


1859 (1941) 


642 (580) 


61.8 


46.3 


(E) 


(0.3) 


(0.3) 


44.2 


35.3 



interactions. From comparisons of these sources and MC 
simulation, the absolute calibration error was estimated 
to be less than ±0.74% (±1.6%) for SK-I (SK-II) period. 

In order to select proton decay signals from the data 
samples, the following selection criteria were applied: (A) 
The number of rings is two or three. (B) One of the rings 
is e-like (/i-like) for p —> e + ir° (p — > /i + 7r°) and all the 
other rings are e-like. (C) For three ring events, 7r° invari- 
ant mass is reconstructed between 85 and 185 MeV/c 2 . 
(D) The number of electrons from muon decay is (1) for 
p — > e + 7r° (p — » /i + 7r°). (E) The reconstructed total mo- 
mentum is less than 250 MeV/c, and the reconstructed 
total invariant mass is between 800 and 1050 MeV/c 2 . 

The detection efficiencies at each selection criterion for 
the p — > e + 7r° and p — > n + ir° search are shown in Tabled] 
The efficiencies were estimated to be 44.6% (43.5%) and 
35.5% (34.7%) for SK-I (SK-II), respectively. The differ- 
ence between the detection efficiencies of SK-I and SK- 
II was only 1%, which implies that SK-II has a similar 
performance with SK-I for these types of proton decay 
events. The inefficiency for proton decay detection is 
mainly due to nuclear interaction effects of pions in 16 O. 
In the simulated proton decay samples, about 37% of 
7r°'s from proton decay in 16 O were absorbed or charge- 
exchanged by interactions with nucleons. Those events 
rarely survive the selection criteria. The efficiency for 
p — > /! + 7T° was lower than that for p — > e + ir° because in 
criterion (D) one electron from muon decay was required 
for the p — ► /i + 7r° search and the detection efficiency for 
electrons from muon decay was approximately 80%. 

The number of remaining events of the atmospheric 
neutrino MC and the data at each selection criterion 
are also shown in Table |H The atmospheric neutrino 
MC, in which neutrino oscillation was accounted for with 
Am 2 = 2.5 x 10 _3 eV 2 and sin 2 26 = 1.0, was normal- 
ized with the observed data using the number of single- 
ring e-like events, which are assumed to have negligi- 
ble neutrino oscillation. Comparison of the number of 
events in TableUdemonstrates that the atmospheric neu- 
trino MC reproduces the observed data well. The num- 
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(stat.)^;f?(sy S .) 



e+7r°. 



ber of background events from atmospheric neutrinos 
in the 91.7+49.2 kiloton-year exposure for p — » e + 7r° 
and p — > (U+7T selection criteria were estimated to be 
0.30 ± 0.04(MC stat.) ± 0.11(sys.) events (= 2.1 ± 0.3 ± 
0.8 events/megaton-year) and 0.34 ± 0.05(MC stat.) ± 
0.12(sys.) events (= 2. 4±0.4±0. 9 events/megaton-year), 
respectively. Approximately 81%(96%) of background 
events for p — > e + 7r° (p — > p + 7r°) were due to charged cur- 
rent neutrino interactions, which consisted of 32%(47%) 
single-pion production, 19%(21%) multi-pion produc- 
tion, and 28%(15%) quasi-elastic scattering (CCQE). In 
most background events from CCQE, a highly energetic 
proton (>1GcV/c) produced by the neutrino interaction 
scatters in water and produces a secondary 7r° which 
makes a similar event signature to proton decay. For the 
systematic error of the background rate, uncertainties of 
event reconstruction, hadron propagation in water, pion- 
nuclear effects, cross sections of neutrino interactions and 
neutrino fluxes were considered. 

The background event rates were also estimated by an- 
other neutrino interaction model, NUANCE [2lj, and the ex- 
perimental data of K2K v beam and the 1-kiloton water 
Cherenkov detector [22|. The expected background rates 
were 1.9±0.7(MC stat.) events/megaton-year by NUANCE 
for both decay modes and 1.63^ 
events/megaton-year by the K2K data for p 
These estimates are consistent with those from our at- 
mospheric neutrino MC. 

We have searched for proton decay signals in fully con- 
tained event samples in SK-I and SK-II. Figure [U(a3) 
and (b3)) and Fig.[2]show distributions of total invariant 
mass and total momentum for the samples satisfying the 
selection criteria except for (E) compared with proton 
decay MC (Fig. [1] (al) and (bl)) and atmospheric neu- 
trino MC (Fig.Q](a2) and (b2)). After applying all the 
proton decay event selection criteria, no candidate events 
for p — > e + 7r° and p — > /i + 7r° were found in the data. 

Because there were no candidates, we calculated lower 
limits on the partial lifetime. In the limit calculation, 
uncertainties on the exposure were negligible. The total 
systematic uncertainty on the detection efficiency was es- 
timated to be 19% for both p — > e + 7r° and p — » /i + 7r° 
and for both SK-I and SK-II. The largest contribution to 
the detection efficiency uncertainty came from uncertain- 
ties in the cross sections for pion- nuclear effects in 16 O; 
these were estimated to be 15% for the p — > e + 7r° and 
p — > p + 7r° detection efficiencies by comparing the pion es- 
cape probability from a nucleus with another model [23j|. 
The systematic error from the uncertainty of the Fermi 
motion was estimated by comparing its momentum dis- 
tribution with the Fermi gas model; this changes the de- 
tection efficiencies by 9%. The fraction of proton decays 
correlated with other nucleons in 16 O also contributes 
an uncertainty that we conservatively set to 100% uncer- 
tainty for a fraction of events; this changes the detection 
efficiencies by 7%. There were no significant differences 
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FIG. 1: Total momentum versus total invariant mass dis- 
tributions for events of proton decay MC (al,bl), 500 year- 
equivalent atmospheric neutrino MC (a2,b2) and data (a3,b3) 
which satisfy the selection criteria of p — > e + ir° except for (E) 
in SK-I(left figures) and SK-II(right figures). The boxes in 
figures indicates the criterion (E). Points in signal boxes of 
the atmospheric neutrino MC (a2,b2) are shown in a larger 
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FIG. 2: Total momentum versus total invariant mass distri- 
butions for events in SK-I(left) and SK-II(right) data, which 
satisfy the selection criteria of p — > /i + 7r° except for (E). The 
box in figures show the criterion (E). 
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between p — > e + 7r° and p — » ^ + 7r° in these uncertain- 
ties. Other sources of uncertainties from energy scale 
stability, particle identification and decay electron find- 
ing efficiency were estimated to be negligible. 

The lifetime limits were calculated by a method based 
on Bayes theorem that incorporates the systematic un- 
certainty These calculations give limits on the partial 



of t/B p 



> 8.2 x 10 33 years 



> 6.6 x 10 33 years at 90% 



lifetime for p — > e 
and p — > /i+7r° of r/B p _ 
confidence level. 

We have reported the results of proton decay 
searches in a 140.9 kiloton-year exposure of the Super- 
Kamiokande detector. The proton decay search perfor- 
mance of SK-II is comparable with that of SK-I. No ev- 
idence for proton decays via the mode p — > e + 7r° and 
p — > /i + 7r° were found and we set limits on the partial 
lifetime of proton. These limits are more stringent com- 
pared with the previous results of 1.6 x 10 33 yearsQ and 
4.7 x 10 32 years 0, respectively. 

We gratefully acknowledge the cooperation of the 
Kamioka Mining and Smelting Company. The Super- 
Kamiokande experiment was built and has been operated 
with funding from the Japanese Ministry of Education, 
Science, Sports and Culture, and the United States De- 
partment of Energy. 
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